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Weyl semimetals are three-dimensional topological states of matter, in a sense that they host 
paired monopoles and antimonopoles of Berry curvature in momentum space, leading to the chiral 
anomaly. The chiral anomaly has long been believed to give a positive magnetoconductivity or 
negative magnetoresistivity in strong and parallel fields. However, several recent experiments on 
both Weyl and Dirac topological semimetals show a negative magnetoconductivity in high fields. 
Here, we study the magnetoconductivity of Weyl and Dirac semimetals in the presence of short- 
range scattering potentials. In a strong magnetic field applied along the direction that connects 
two Weyl nodes, we find that the conductivity along the field direction is determined by the Fermi 
velocity, instead of by the Landau degeneracy. We identify three scenarios in which the high-held 
magnetoconductivity is negative. Our hndings show that the high-held positive magnetoconductivity 
may not be a compelling signature of the chiral anomaly and will be helpful for interpreting the 
inconsistency in the recent experiments and earlier theories. 

PACS numbers: 75.47.-m, 03.65.Vf, 71.90.-|-q, 73.43.-f 


I. INTRODUCTION 


Topological semimetals are three-dimensional topolog¬ 
ical states of matter. Their band structures look like 
three-dimensional analogue of graphene, in which the 
conduction and valence energy bands with linear disper¬ 
sions touch at a finite number of points, i.e., Weyl nodes 
[1]. The nodes always occur in pairs and carry oppo¬ 
site chirality. One of the topological aspects of Weyl 
semimetals is that they host pairs of monopole and anti¬ 
monopole of Berry curvature in momentum space [2, 3] 
(see Fig. 1), and the fluxes of Berry curvature flow from 
one monopole to the other. In the presence of both a 
magnetic field and an electric field along the direction 
that connects two monopoles, electrons can be pumped 
from one monopole to the other, leading to the Adler- 
Bell-Jackiw chiral anomaly [4-6] (also known as trian¬ 
gle anomaly). Recently, angle-resolved photoemission 
spectroscopy (ARPES) has identified the Dirac nodes [7] 
(doubly-degenerate Weyl nodes) in (Bii_ 2 ;In 2 ,) 2 Se 3 [8, 9], 
NasBi [10-13], and Cd 3 As 2 [12, 14-17] and Weyl nodes in 
TaAs [18-21]. Also, scanning tunneling microscopy has 
observed the Landau quantization in Cd 3 As 2 [22] and 
TlBiSSe [23]. 

While the chiral anomaly is well established in momen¬ 
tum space, it becomes a challenging issue how to detect 
the effects of the chiral anomaly or relevant physical con¬ 
sequences. This has been attracting a lot of theoretical 
efforts, such as the prediction of negative parallel mag¬ 
netoresistance [24-27], proposal of non-local transport 
[28], using electronic circuits [29], plasmon mode [30], 
etc. In particular, whether or not the chiral anomaly 
could produce measurable magnetoconductivity is one of 
the focuses in recent efforts. This has inspired a num¬ 
ber of transport experiments in topological semimetals 
Cd 3 As 2 [31-36], ZrTes [37], NbP [38], NagBi [39], and 
TaAs [40, 41]. The chiral anomaly has been claimed 
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FIG. 1. Left: the energy spectrum of a Weyl semimetal as 
a function of at k,,, = ky — 0. Right: the vector plot of 
the Berry curvature of a Weyl semimetal, {k^, ky,kz) is the 
wave vector. The analytic expression of the three-dimensional 
Berry curvature is given in Eq. (3). The parameters are 
Mo = Ml = A = 1, so fcc = 1. 


to be verified in several different topological semimet¬ 
als, including BiSb alloy [42], ZrTes [37], TaAs [40, 41], 
and Na 3 Bi [39], in which similar magnetoconductivity 
behaviors are observed when the magnetic field is ap¬ 
plied along the conductivity measurement direction: (i) 
In weak fields, a —'/B negative magneto conductivity is 
observed at low temperatures, consistent with the quan¬ 
tum transport theory of the weak antilocalization of Weyl 
or Dirac fermions in three dimensions [42, 43]. (ii) In in¬ 
termediate fields, a positive magnetoconductivity is 
observed, as expected by the theory of the semiclassical 
conductivity arising from the chiral anomaly [26, 27, 44]. 
(iii) In high fields, the magnetoconductivity is always 
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negative in the experiments. However, in the strong-held 
limit, a positive magnetoconductivity is expected in ex¬ 
isting theories, also as one of the signatures of the chiral 
anomaly [24-26, 45]. 

In this work, we focus on the high-held limit, and 
present a systematic calculation on the conductivity of 
topological semimetals. Beyond the previous treatments, 
we start with a two-node model that describes a pair of 
Weyl nodes with a hnite distance in momentum space. 
Moreover, we fully consider the magnetic held depen¬ 
dence of the scattering time for electrons on each Landau 
level, and obtain a conductivity formula. The efforts lead 
to qualitatively distinct results compared to all the pre¬ 
vious theories. We hnd that the conductivity does not 
grow with the Landau degeneracy but mainly depends 
on the Fermi velocity. The magnetoconductivity arises 
from the held dependence of the Fermi velocity when the 
chemical potential is tuned by the magnetic held. Based 
on this formula and the model, we hnd that although 
the positive magnetoconductivity is also possible, three 
cases can be identihed in which the magnetoconductivity 
is negative, possibly applicable to those observed in the 
experiments in high magnetic helds. 

The paper is organized as follows. In Sec. II, we in¬ 
troduce the two-node model and show how it carries all 
the topological properties of a topological semimetal. In 
Sec. Ill, we present the solutions of the Landau bands 
of the semimetal in a magnetic held applied along the z 
direction (the two Weyl nodes are separated along this 
direction). In Sec. IV, we calculate the z-direction mag¬ 
neto conductivity in the presence of the short-range delta 
scattering potential. In Sec. V, we discuss various sce¬ 
narios that the negative or positive magnetoconductivity 
may occur. In Sec. VI, we present the transport in the 
X — y plane, including the x-direction conductivity and 
the Hall conductance. Finally, a summary of three sce¬ 
narios of the negative magnetoconductivity is given in 
Sec. VH. The details of the calculations are provided in 
Appendices A-D. 


to two separate local models 


H± = M± • cr. 


( 2 ) 


H± = Mj- • CTwith Mj- = (^Akx, Aky,^2Mikc.k^ and 

{kx,ky,kz) the effective wave vector measured from the 
Weyl nodes. 

The topological properties in H can be seen from the 
Berry curvature [47], H(k) = Vk x A(k), where the Berry 
connection is dehned as A(k) = i (u(k)j Vk [^(k)). For 
example, for the energy eigenstates for the + band lu(k)) 
= [cos(0/2),sin(0/2)e®‘^], where cos0 = Mk/E+ and 
iaxiip = ky/kx- The three-dimensional Berry curvature 
for the two-node model can be expressed as 


n(k) = 


A^Mi 

~W 


kzkx,kzky, ^{k^ k^ kx ky) 


( 3 ) 


When MqMi > 0, there exist a pair of singularities 
at (0,0,±A:c) as shown in Fig. 1. The chirality of a 
Weyl node can be found as an integral over the Fermi 
surface enclosing one Weyl node (l/27r)^ H(k) • (iS(k), 
which yields opposite topological charges =Fsgn(Mi) at 
Eke, corresponding to a pair of “magnetic monopole 
and antimonopole” in momentum space. For a given 
kz, a Chern number can be well defined as Ueikz) = 
— (I/27r) ff dkxdkyfi(k) ■ z to characterize the topological 
property in the kx-ky plane, and ndkz) = —^[sgn(Mo — 
Mikl) + sgn(Mi)] [48]. For MqMi > 0, ndkz) = 
—sgn(Mi) for —kc < kz < kc, and ndkz) = 0 for other 
cases [49] . The nonzero Chern number corresponds to the 
fc^-dependent edge states (known as the Fermi arc) ac¬ 
cording to the bulk-boundary correspondence [50] . Thus 
the two-node model in Eq. (1) provides a generic descrip¬ 
tion for Weyl semimetals, including the band touching, 
opposite chirality, monopoles of Berry curvature, topo¬ 
logical charges, and Fermi arc. In the following we shall 
focus on the topological case of MqMi > 0. 


II. MODEL AND ITS TOPOLOGY 

A minimal model for a Weyl semimetal is 

H = A{kxCfx + kyay) + M\,,az, (I) 

where a are the Pauli matrices, Aik = Mq — Mi{kx + 
ky + k1), {kx,ky,kz) is the wave vector, and A, Mq/i 
are model parameters. This minimal model gives a 
global description of a pair of Weyl nodes of opposite 
chirality and all the topological properties. It has an 
identical structure as that for A-phase of ^He super¬ 
fluids [46] The dispersions of two energy bands of this 

model are E± = iy^AI^ + A'^{kx + k^), which reduce 
to E± = E\Mq — Mik'd at kx = ky = 0. If MqMi > 0, 
the two bands intersect at (0, 0, Ekd) with kc = i/MqJMi 
(see Fig. 1). Around the two nodes (0, 0, Eke), H reduces 


III. LANDAU BANDS 

In a magnetic field along the z direction, the energy 
spectrum is quantized into a set of ID Landau bands dis¬ 
persing with kz [see Fig. 2 (a)]. We consider a magnetic 
field applied along the z direction, B = (0,0,5), and 
choose the Landau gauge in which the vector potential 
is A = (—j/5,0, 0). Under the Pierls replacement, the 
wave vector in the Hamiltonian in Eq. (1) is replaced by 
the operator 

pR 

^^{kx - ^y,-idy,kz) (4) 

kx and kz are still the good quantum numbers as the in¬ 
troduction of the gauge field does not break the transla¬ 
tional symmetry along the x and z direction. Introducing 
the ladder operators [51, 52], k'^ E k'^ w(a'^a -I- 1/2), 
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FIG. 2. The energies of Landau bands of the minimal global 
model for Weyl and Dirac semimetals in a magnetic field B 
applied along the 2 direction, as functions of the wave vector 
kz- The parameters; Mq = 0.05 eV, Mi = 5 eV-nm^, A = 1 
eV-nm, and B — 1 Tesla. The Zeeman energy is not included. 


k+ —)■ (•\/2/€B)a'^, k- —)■ {\/2liB)a, where the mag¬ 
netic length £b = i/h/eB and the ladder operators 
a = -[(y - i%ka:)/(.B + £Bdy\l\/2 and a'^ = -[{y - 
£‘Bkx)/£B — £Bdy]lV2 [51, 53], then we can write the 
Hamiltonian in terms of the ladder operators. 


H{k) = 


Mfc Ak- 


' Ma 

r]a 

1 

1 

+ 

_1 


rja^ 

-Ma_ 


( 5 ) 


where Aia = Mq — Mik^ — a;(a^o -I- 1/2), uj = 
and rj = iflAjiB- With the trial wave functions 
(cijr^— 1), C 2 |j^))^ for = 1,2,... (later denoted as > 1) 
and (0, |0))^ for 1 ^ = 0, where ly indexes the Hermite 
polynomials, the eigen energies E can be found from the 
secular equation 


det 


A4i,+ujj2 — E 

r)i/v —M^+uj/2 — E 


( 6 ) 


for V > and —M-i, -\- oj/2 — E = 0 toT 1 / = 0, where 
Ml, = Mq — Mikl — uv. The eigen energies are 

E-'^f = w/2 ± + V > 1 

El^ = uj/2 - Mo -b Mikl, u = Q. (7) 

This represents a set of Landau energy bands (ly as band 
index) dispersing with k^, as shown in Fig. 2. The eigen 
states for > 1 are 


I^ 1, kx , ^2 , “b) — 


cos — 1) 


1 ^ ^ 1 , kxj kzi ) — 

and for ;/ = 0 is 

\v = 0,kx,kz) = 


sm 


k-1) 


I kx ikz') 1 

\kx,kz), 


0 

| 0 ) 


\kx,kz), 


( 8 ) 


( 9 ) 


where cos 9 = M^l\/Ml + vrf^ and the wave functions 
'0j/,fc,.fe,(r) = (rk, kx-, kz) are found as 


(l*) - 


Ciz 




(y-ypr 


VLxLz£b 


HA 


y-yo. 




( 10 ) 


where Cj, = 1/i/iA^ a/tt, LxLz is area of sample, the 
guiding center yo = kxt%, Hi, are the Hermite polynomi¬ 
als. As the dispersions are not explicit functions of kx, 
the number of different kx represents the Landau degen¬ 
eracy Nl = \/2'k(?b = eB/h in a unit area in the x-y 
plane. 

This set of analytical solutions provides us a good 
base to study the transport properties of Weyl fermions. 
In the following, we will focus on the quantum limit, 
i.e., only the ly = 0 band is on the Fermi surface [see 
Fig. 4 (b)]. 


IV. Z-DIRECTION SEMICLASSICAL 
CONDUCTIVITY 

A. Argument of positive magnetoconductivity 

When the Fermi energy is located between the two 
states of \iy = l,kx,kz,E), all the bands for |i/ > 
1, kx, kz, -b) are empty and all the bands \iy > 1, kx, kz, —) 
are fully occupied. Only the band of = 0 is partially 
filled. In this case the transport properties of the sys¬ 
tem are dominantly determined by the highly degener¬ 
ate V = Q Landau bands [the red curve in Fig. 2 (a)] 
. It is reasonable to regard them as a bundle of one¬ 
dimensional chains. Combining the Landau degeneracy 
Nb, the z-direction conductance is approximately given 
by 


CTzz = NlCJib, (11) 


where ctid is the conductance for each one-dimensional 
Landau band. 

If we ignore the scattering between the states in the de¬ 
generate Landau bands, according to the transport the¬ 
ory, the ballistic conductance of a one-dimensional chain 
in the clean limit is given by 

c^iD = (12) 


then the conductivity is found as 


eB 


(13) 


which is is linear in magnetic held B, giving a positive 
magnet o conductivity. 

In most measurements, the sample size is much larger 
than the mean free path, then the scattering between the 
states in the Landau bands is inevitable, and we have to 
consider the other limit, i.e., the diffusive limit. Usually, 
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the scattering is characterized by a momentum relaxation 
time T. According to the Einstein relation, the conduc¬ 
tivity of each Landau band in the diffusive limit is 

CTiD = e^A^iDUpT, (14) 


where vp the Fermi velocity and the density of states for 
each ID Landau band is A^id = l/T^hvp, then 


eBvpT 
h ph 


(15) 


If and r are constant, one readily concludes that the 
magneto conductivity is positive and linear in B. 

According to Nielsen and Ninomiya [24], to illustrate 
the physical picture of the chiral anomaly, they started 
with a one-dimensional model in which two chiral en¬ 
ergy bands have linear dispersions and opposite veloci¬ 
ties. An external electric field can accelerate electrons 
in one band to higher energy levels, in this way, charges 
are “created”. In contrast, in the other band, which has 
the opposite velocity, charges are annihilated. The chiral 
charge, defined as the difference between the charges in 
the two bands, therefore is not conserved in the electric 
field. This is literally the chiral anomaly. As one of the 
possible realizations of the one-dimensional chiral system, 
they then proposed to use the v = 0 Landau bands of a 
three-dimensional semimetal, and expected “the longi¬ 
tudinal magneto-conduction becomes extremely strong”. 
In other words, the magnetoresistance of the 0th Landau 
bands in semimetals is the first physical quantity that was 
proposed as one of the signatures of the chiral anomaly. 

Recently, several theoretical works have formulated the 
negative magnetoresistance or positive magnetoconduc¬ 
tivity in the quantum limit as one of the signatures of 
the chiral anomaly [26, 45], much similar to those in Eqs. 
(13) and (15). In both cases, the positive magnetocon¬ 
ductivity arises because the Landau degeneracy increases 
linearly with B. However, in the following, we will show 
that if uf and r also depend on the magnetic field, the 
conclusion has to be reexamined. 


B. Green function calculation 


where —e is the electron charge, V = L^LyL^ is the 
volume with the length along the x direction and so 
on, Ug = dE^^/hdkz = 2Mik^/h is the velocity along 
the z direction for a state with wave vector kz in the 
p = 0 band, = l/{Ep — E^^ ± ihj2T°) is the re¬ 

tarded/advanced Green’s function, with the lifetime 
of a state in the p = 0 band with wave vector kx and 
kz- Usually, in the diffusive regime, one can replace 
by ^T°J{Ep - E°J. ^ However, in one di¬ 
mension, to correct the van Hove singularity at the band 
edge, we introduce an extra correction factor A, so that 
GokG^k^ = -E^J- As shown in Appendix 

A, A —>■ I (0) if the Fermi energy is far away from (ap¬ 
proaching) the band edge. Now the conductivity formula 
can be written as 

^-.0 = ^ E {vk.r^ErlSiEp - El). (17) 

kz ,kx 


The delta function S{Ep — El) = ''^bere Vp = 

2\Mi\kp/h is the absolute value of the Fermi velocity of 
the p = 0 band with kp the Fermi wave vector. This 
allows us to perform the summation over kz, then 


(T 


SC _ 

zz,0 


VAr° 

hLxLy^ 


(18) 


The summation over kx is limited by the Landau degen¬ 
eracy, finally we can reduce the conductivity formula to 


cr 


SC _ 

22,0 ~ 



h PLy 



dkxT^E. 


(19) 


The scattering time depends on the wave packet of 
the Landau levels in band 0 and is a function of magnetic 
field. It can be found from the iteration equation under 
the self-consistent Born approximation (see Appendix B 
for details) 


= E 




0,0 

kxkz.k'kL 


'-)W2rl 


KM. 


■, [{Ep-El)^K{h/2r^)^\' 


( 20 ) 


Now we are ready to present the conductivity in the 
presence of the magnetic field when the Fermi energy is 
located near the Weyl nodes. The temperature is as¬ 
sumed to be much lower than the gap between bands 
H- and I—, i.e, kpT 2-\/2A/£b. In this case all the 
Landau levels of are fully occupied while the p — 0 
band [the red curve in Fig. 2 (a)] is partially filled. Since 
El is only a function of kz, and independent of kx, the 
system can be regarded as a bundle of highly degener¬ 
ate one-dimensional chains. Along the z direction, the 
semiclassical Drude conductivity can be found from the 
formula [54] 

'"-.0 = ^ E ( 16 ) 

kz ikx 


where Ul\ j,, represents the scattering matrix ele¬ 
ments, (...) means the average over impurity configura¬ 
tions. The conductivity in semimetals in vanishing mag¬ 
netic field has been discussed within the Born approxi¬ 
mation [55]. 

In this work, we consider only the short-range delta 
scattering potential. The delta potential takes the form 

U(r)=^u,<5(r-R,), (21) 

i 

where Ui measures the strength of scattering for an 
impurity at R^, and the potential is delta correlated 
{U {v)U {r')) = Uimp<5(r — r'), where Uimp is a field- 
independent parameter that is proportional to the im¬ 
purity density and averaged field-independent scattering 
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strength. Using the wave function of the v = 0 band, we 
find that (see Appendix C) 


r0,0 


^ I kz ,k' k' 


I2\ _ Virnpexp[-e%{kj; - K)^/2] 


1 = 


^ z^B '\/27r 


( 22 ) 


and in the strong-field limit {B —^ oo), 


^A = 






imp 


(23) 



which gives the conductivity in the strong-held limit as 


(hvp)'^ 
^ l^mp 


(24) 


Notice that the Landau degeneracy in the scattering time 
cancels with that in Eq. (19), thus the magnetic held de¬ 
pendence of g is given by the Fermi velocity Vp. This 
is one of the main results in this paper. When ignor¬ 
ing the magnetic held dependence of the Fermi velocity, 
a B-independent conductivity was concluded, which is 
consistent with the previous work in which the velocity 
is constant [25]. Later, we will see the magnetic held 
dependence of the Fermi velocity can lead to different 
scenarios of positive and negative magnetoconductivity. 


V. SCENARIOS OF NEGATIVE AND POSITIVE 
Z-DIRECTION MAGNETOCONDUCTIVITY 


A. Weyl semimetal with fixed carrier density 


FIG. 3. (a) The z-direction conductivity of the 0-th Landau 
band of the Weyl semimetal as a function of the 2 : direction 
magnetic field B, for different values of the carrier density no 
(in units of 10^^/cm^). The lower bound of the field for each 
curve is determined by the bottom of band U- for no > 0 
or the bottom of band 0 for no < 0. (b) The characteristic 
Be defined in Eq. (27) as a function of no. Parameters: 
Mo = 0.05 eV, Ml = 5 eV-nm^, A = 1 eV-nm, Vimp = 10 
(eV)^-nm^. Using no = {Stt/ 3)E%/2MikcA'^{2 tt)^, a carrier 
density of no x 10^^/cm®corresponds to a Fermi energy of 
about Ef =144x®y^ho meV. 


Fig. 3 (b)]. g is constant for the undoped case of 
rig = 0, and 


h Ump 


(28) 


is the conductivity of the undoped case, and is indepen¬ 
dent of magnetic held. Thus the magneto conductivity is 
always negative in the electron-doped case while always 
positive in the hole-doped regime as shown in Fig. 3 (a). 


In a strong held the Fermi velocity or the Fermi energy 
is given by the density of charge carriers and the magnetic 
held [56]. We assume that an ideal Weyl semimetal is 
the case that the Fermi energy crosses the Weyl nodes, 
all negative bands are fully hlled and the positive bands 
are empty. In this case hv^p = 2Mikc- An extra doping 
of charge carriers will cause a change of electron density 
no(> 0) in the electron-doped case or hole density no(< 
0) in the hole-doped case. The relation between the Fermi 
wave vector and the density of charge carriers is given by 

ng = 2Np X (25) 

ZTT 

This means that the Fermi wave vector is determined by 
the density of charge carriers rig and magnetic held B, 

k^p = kc + 'KUph/eB (26) 

or kp = kc + 27r^ng£^. Thus the Fermi velocity is also a 
function of B, Hvp = 2Mikp, and 


B. Weyl semimetal with fixed Fermi energy 

In the case that the Fermi energy is hxed, (hvp)'^ = 
AMi{Ep — eMiB/h + Mg), and we have 

e2 4M^{Ep - eM,B/H + Mg) 

^Z2,0 7 TA I (29) 

7 imp 

then the magnetoconductivity is always negative and lin¬ 
ear in B. 


C. Dirac semimetal 

If the system has time-reversal symmetry, we may have 
a Dirac semimetal, instead of Weyl semimetal, and all 
Weyl nodes turn to doubly-degenerate Dirac nodes. A 
model for Dirac semimetal can be constructed [12] by 
adding a time-reversal partner to Eq. (1) 


^zz,0 - 


1 -b sgn(ng) 


Be 

B 


n 2 


(27) 


L^Dirac — 


iL(k) 0 
0 M*(-k) 



' A, 

0 1 

+ CTz (g) 

0 

-1 

P, 

<1 


■ (30) 


where the characteristic held Be = TT jngj/i/eAic. A typi- In the second term, the 2 -direction Zeeman energy 
cal order of B^ is about 10 Tesla for ng of 10^^/cm^ [see ^s/p = 9s/pk‘BB/2 is also included, where is the 
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g-factor for the s/p orbital [22] and ps is the Bohr 
magneton. Fig. 2 (b) shows the Landau bands of both 
iL(k) and iL*(—k) in the ^-direction magnetic field. The 
Landau bands of the Dirac semimetal can be found in 
a similar way as that in Sec. III. Now there are two 
branches oi v = 0 bands, with the energy dispersions 
= w/2 + Ap — Mq + Mik'^ and = — uj/2 — + 

Mq — Mikl for iL(k) and iJ*(—k), respectively. They 
intersect at k^ = ±-\/[Mo — (w + A^ + Ap)/2]/Mi and 
energy (Ap — As)/2, and with opposite Fermi velocities 
near the points. In the absence of inter-block velocity, 
the longitudinal conductance along the z direction is ap¬ 
proximately a summation of those for two independent 
Weyl semimetals. 

First, we consider the Fermi energy cross both bands 
0 t and 0 1. Using Eq. (29), the z-direction conductivity 
is found as 




FIG. 4. (a) The i-direction conductivity of the 0-th Landau 
band of the Weyl semimetal as a function of the z direction 
magnetic field B, for different values of the carrier density 
no (in units of 10^^/cm^). (b) Schematic of the second-order 
processes that contribute to the x direction conductivity. Pa¬ 
rameters: Mo = 0.05 eV, Ml = 5 eV-nm^, A = 1 eV-nm, 
Uimp = 10 (eV)^-nm®. 


cr 


SC 

zz,0 


= a 


sc 

22,Ot 


I _SC 

^zzM 


eMiB 

-^[Mo-^ 

^ Mmp ^ 


f^B{9p + gs)B ^ ( 31 ) 


VI. LONGITUDINAL AND HALL 
CONDUCTIVITIES IN X-Y PLANE 


A. ^-direction conductivity 


or using defined in Eq. (28), 


a 


SC _ 

zz,0 ~ 


2crAr[l 


eB 

hk/ 


k-B{gp + gs)B 

4Mo 


(32) 


In this case we have a negative linear B magnetocon¬ 
ductivity, when the Fermi energy crosses both and 

eZ 

With increasing magnetic field, the 0 f bands will shift 
upwards and the 0 bands will shift downwards. Be¬ 
yond a critical field, the Fermi energy will fall into either 
0 t or 0 bands, depending on whether the carriers are 
electron-type or hole-type. If the carrier density is fixed, 
the Fermi wave vector in this case does not depend on k^ 
as that in Eq. (26), but 


In the X — y plane normal to the magnetic field, the 
longitudinal conductivity along either the x or y direction 
is negligibly small as the effective velocity 


V 


X 

0 


dEt 

hdkx 


(35) 


Nevertheless, a non-zero longitudinal conductivity along 
the X direction can be found as 

= ^ E Re(GoXmG^t<o). (36) 

kx,kz 


where the Green’s functions of band 1± are = 

l/iEp — E]/^ ± ih/2T//^f,^), and the inter-band velocity 


0 _ Trno/i 


(33) 


or kPp = 27r^no^5. case, with increasing magnetic 

field, the Fermi energy will approach the band edge and 
the Fermi velocity always decreases. Using Eq. (24), 


cr 


SC _ 

22,0 ~ 


~h 


(34) 


which also gives negative magnetoconductivity that is 
independent on the type of carriers. 

Note that in the Weyl semimetal TaAs with broken in¬ 
version symmetry, where the Weyl nodes always come in 
even pairs because of time-reversal symmetry [18-21], the 
situation is more similar to that for the Dirac semimetal 
and the magnetoconductivity does not depend on the 
type of carriers and may be described by a generalized 
version of Eqs. (32) and (34). 


Ki+ = 
K.i- = 


Sin - 


V2Mi . e 

^B 

y/2Mi 


A cos 


^B 


COS —^ -I- A sm —^. 


(37) 


Note that for the Landau bands generated by the z- 
direction magnetic field, the leading-order a:-direction ve¬ 
locity is the inter-band velocity v/j that couples band 
0 with bands 1±, and the scattering times of 

band 1± are due to virtual scattering going back and 
forth between band 0 [see Fig. 3 (d)], so cr^x,o indeed 
stems from second-order processes and therefore is much 
smaller than crf^ q that arises from first-order processes. 
We find that cr®'), q = , where (see Appendix 

D for details) 


,o.i± 


e h 


,A6{Ef-E°) h 


V 2(Ef - E/^y 


.(38) 


kx ,kz 


k:xkz 
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At this stage, we have the same form as Eq. (28) in 
the paper by Abrikosov [56]. If the Hamiltonian is re¬ 
placed hy H = vk-a and the scattering time is evaluated 
for screened charge impurities under the random phase 
approximation, an a;-direction 1/B magnetoconductivity 
can be found, leading to the quantum linear magnetore¬ 
sistance. In the present case, the scattering time is found 
as 


In the state oi v = 0, the energy dispersion is corrected 
to — eEyi'^kx as {y) = l^k^- This energy correction 
leads to a velocity shift along the x direction, 

_ 1 diEl - eEyilk,) eEyll 

"■ = -h -at-= 

For each kz , this Vx leads to a quantized Hall conductance 


= (=>9) 

Note that is proportional to so the Landau de¬ 
generacy l/27rl!^ from the conductivity formula cancels 
with that from the scattering time. However, q is 
inversely proportional to the scattering time then 

the effect of the Landau degeneracy actually is doubled, 
and finally we arrive at 


A^ H,i+ikz)\^ . 2^*^ 

(40) 

where one replaces sin by cos for . In both the 
electron- and hole-doped regimes, the magnetoconduc¬ 
tivity is always positive as shown in Fig. 3 (c). 

In the strong-held limit, sin^(0^^/2) —)• I if Mi > 0, 
kiVQ^y,+ —>■ Mi\/2/(.b according to Eq. (37), and {Ep — 
E],^)'^ = AMi/£g according to Eq. (7), then 


,0.1-K _ 




h 2TT^i% 


^ 1 _ 

Ey Ey 27r^^ h 


(46) 


The total Hall conductance is found by integrating over 
kz up to the Fermi wave vector k^p, and 


<^yx — ( 47 ) 

In particular, kp = kc + irnoh/eB for Weyl semimetals 
with a hxed carrier density no, and a Hall conductance 
is found as 


^yx 


enp 

B 



(48) 


The hrst term is attributed to the classical Hall ef¬ 
fects, and the second term comes from the non-zero 
Chern number of the fully hlled low energy bands of 
—kc < kz < kc- This is consistent with the calculation 
by using the Kubo formula for the Hall conductivity [60] . 


a 


0,1 

XX 


u A 

^ ^ imp / \2 

h hvp 


(41) 


VII. SUMMARY 


(1) In Weyl semimetals with a hxed carrier density, 
the magnetic held will push the Fermi wave vector to 
^F = near which A —>■ 1 and hvp = 2Mikc{l + 
sgn(no)Hc/H), then 


0,1 _ ^ _ eVirapB _ 

h IGn'^hMikcil + sgn{no)Bc/B)'^ 

In the limit that B ^ Be, cr^’x increases linearly with B. 

(2) In Dirac semimetals, the magnetic held pushes the 
Fermi energy to the band edge, using A/z;^ in Eq. (A8), 


^0,1 


\ ^ ^1/3 

h 47 t2£% [aMiVi^J 


(43) 


B. Hall conductivity 

In the presence of the z-direction magnetic held, a Hall 
conductance in the x — y plane can also be generated 
[49, 57-59]. The correction of an electric held Ey to the 
model Hamiltonian is the potential energy, 


(44) 


We present a conductivity formula for the lowest Lan¬ 
dau band in a semimetal in the presence of short-range 
delta scattering potentials. The conductivity depends on 
the square of the Fermi velocity, instead of the Landau 
degeneracy. Based on this mechanism and the model that 
describes two Weyl nodes with a hnite spacing in momen¬ 
tum space, we hnd three cases that could give a negative 
magneto conductivity in the strong-held limit, (i) A Weyl 
semimetal with a hxed density of electron-type charge 
carries [Eq. (27)]. (ii) A Weyl semimetal with a hxed 
Fermi energy [Eq. (29)]. (iii) A Dirac semimetal or a 
Weyl semimetal with time-reversed pairs of Weyl nodes 
[Eqs. (32) and (34)], with a 1/H^ dependence. These 
formulas are valid as long as the Fermi energy crosses 
the V = 0 Landau bands. Our theory can be applied to 
account for the negative magneto conductivity observed 
experimentally in various topological semimetals in high 
magnetic helds, such as BiSb alloy [42], TaAs [40, 41], 
and NaaBi [39]. In this way, we conclude that a positive 
magneto conductivity (or negative magnetoresistance) in 
the strong-held limit is not a compelling signature of the 
chiral anomaly in topological semimetals. Our theory can 
also be generalized to understand the magneto conductiv¬ 
ity in non-topological three-dimensional semimetals. 


AV = -eEyy. 
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Appendix A: About the correction factor A 

For the v = Q band, we need to deal with the imaginary 
part of the Green’s function 


h 


[M,kl - {Ef + Mo- cc/2)]2 + (^)2 ’ 


(Al) 


In this work, we assume Mq, Mi, E f,uj > 0- In th is 
case, we can write a = y/Mi, b = \/Ef + Mq — uil2 = 
Hvp/2-/M[, and c = A widely used approxima¬ 

tion is that 


dx 


l/c 


(a2a;2 — 52)2 -p I/c2 


dxTrS{a^x^ — b^). 

— 00 

(A2) 

However, this leads to unphysical van Hove singularities 
at the band edges. We correct this approximation with 
an extra factor A, so that 


dx 


1/c 


(a 2 a ;2 — 52)2 -(- l/c 2 


/ OO 

dxTTAS{a^x^ — b^). 

-00 

(A3) 


The form of A can be found as follows. First, the integral 
can be found as 


dx 


lA 


Ac ^ Ag2 + 


{a^x'^ — b'^)'^ + l/c^ a y/2 vAHAT 


(A4) 

where = b^c. On the other hand, using the property 
of the delta function 

/ OO 

dxd^a^x^ — b^) = irA/ab. (A5) 

-OO 


So 


A = 


(A6) 


In the limit C 1, A —)■ 1, while in the other limit 

C< 1, 


A 2 {fiVp)^T° 

lim A^ = ^ 


(A7) 


c<i 4Mi h ’ 

so this correction is necessary near the band bottom, 
where Ep + Mq — ojf 2 = 0, and can be ignored when 
the Fermi energy is far away from the band bottom. 

Now we evaluate A/Hvp in the band bottom limit. 
Combining Eqs. (23) and (A7), we arrive at 

hm ^ ' oc (A8) 

C«i V4AflUmp/ 


In the self-consistent Born approximation, the full 
Green function is written as 




1 


Ef-EI,-ER{Ef) 


(HI) 


where the self-energy 

E^k.kSEp) = Y, {\UlY^,^,f)^^KK^EF). (B2) 

K.K 

The real part of the self-energy can be absorbed into the 
definition of the chemical potential, we only need the 
imaginary part, 


1 = 




E 


le- 




LMbV^ (Ef - Elf +12 ■ 


(B3) 


where we have used Eq. (22) and suppressed the A 

and kz dependence of X in the strong-held limit. Us- 

; „ JL ■'P _ r°° 1 _ A»/2^b _ 

2^k'^ ~ J -00 27V ’ “ i-Lyj2l\ 27v ~ 

h Ienkl-Lj2ts ^ - K), and in the 

strong-held limit, ± —>■ ± 00 , 


Ump 


dki 


2MIJ_^ 2f {Ep-Elf+X^- 
Using E^^ = uj/2 — Mq + Mf, 

r dk' X 


(B4) 


X = 


2fI% 2f [Mi{kf -Ef+ a;/2 - Mo]2 + 12 ' 

(B5) 

We consider Ep + Mq > a;/2, so the integral can be 
written into 


1 




imp 


dx 


lA 


c (27r)2^^ (a2a;2 — 52)2 + l/c2 


, (B6) 


with a = y/Mi, b = tJEf + Mq — w/2, c = IjX. Using 
Eq. (A5), 


f 


dx 


1/c 


A 

= ^A, 


(a2a;2 — bf + l/c^ ab 


then 




2r° c (27r)2£^ ab 


^ A - ^ 


27r.^^ hvl 


(B7) 


(B8) 


Appendix C: Scattering matrix elements 

We calculate a general form of the scattering matrix 
element where 

= J dr'{a,k'z„kf)U{r'){r'\P,kz:,kz), 

^fcE,/ci/c' = J dr{'y,ka„kz\r)U{r){r\S,k'y„k'^). (Cl) 
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Using the wave functions in Eq. (10) and the identity 

(/ rfr / dr'f{r)f{r')U{r)U{r')) 

= JdrJ dr'f{r)fir'){U{r)U{r')), (C2) 

and {U{r)U{r')) = Uimp<5(r — r'), we have 


/rra,/3 . 


(v-va)^ (y-yp)^ 


VimpCgCfsC^Cs f 

Limi J 

IB ^B ^B 


(C3) 


After performing the integral, we have 


/7-r“./3 \ — TaP~iS (nA^ 


where we have defined a dimensionless integral 


2 

= CgCpC^Cse-^ / dXe-^^"Hg{X + -) 

J — oo ^ 

y,np{x-^)n^{x-^)Hs{x + ^), (C5) 


and K = {yo — y'p)/^B- Using Eqs. (C4) and (C5), it is 
straightforward to find Eq. (22). 


Appendix D: Calculation of ^-direction conductivity 


Now we evaluate 


(J^d± — 


e^h 


and the velocities along x direction are found as 

,Miy/2 . 91 oi 

2 ' 2 


Swo.p+ = sin ^ + Acos -^)(5^.i 






= (- 


Miy/2 


et et 

-cos ^ + Asin 

B 2, I 


(D3) 


then 


, 0.1 + 


27rU 


E2|' 


h 


'0,1+1 


kx^kz 


[{EF-Ely + {^r] 


2t; 


[(EB-El+r + i^y] 


(D4) 


We assume that the Fermi energy crosses only the z/ = 0 
Landau band, in this case, 


h 

2t0 


[[Eb-eIY + {^Y] 


7tAS{Ef-EI). (D5) 


Note that different from Abrikosov’s approximation in 
Eqs. (27)-(28) of Ref. 56, an extra correction factor A is 
added. Then 


(J. 


0 . 1 + ^ 

V ^ 


AS(Eb-EI)^ 


kx,kz 


[Ef-eI+y + {^) 


(D6) 


Using i ^ S{Ef - 

lO.* 


E 


where kz = kp are the roots of 


Lz ^kz J 27r ’ L, 

0 N _ <^(fcf'‘-fc 

kj 

Ep = and for the p = 0 band, = —kp~ = kp, 
and considering \Ep-El+\:s> +++ 


,0.1 + 


e 

h 


I^<i+(fc0l' 


fcifcz 

1 


\hv°p\ [{Ep-El+)^] 


k —k^ 


1 y.h 


yy J-Ly/2el 


27r 


A- 


. 1 + 

‘ kxkz 


(D7) 


k^=k° 


This is evaluated in the numerical calculation. Because 
the 1+ band does not cross the Fermi energy, Tj 


1 + 
fc, k. 


IS 


27rU ^ 

kx ,k. 


2Re q) , (Dl) mainly contributed by the virtual scattering processes 


with the p = 0 band, i.e., t}\ —)• and 


where the Green’s functions 



n 


^l+e^O 


'fe,/Cz 


= 271 ^ {\UlX,k',^k'f)ESiEp - El,). (D8) 


K,K 


Using Eq. (C4) and j)i=i = 2/^+, we 

then find Eqs. (39) and (40). Similarly, in the strong- 


(D2) field limit, one can obtain Tj 


I and cr°d , 

Z Uz 
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